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ABSTRACT 

Context. Thanks to remarkable progress, radial velocity surveys are now able to detect terrestrial planets at habitable 
distance from low-mass stars. Recently, two planets with minimum masses below 10 Mohave been reported in a triple 
system around the M-type star Gliese 581. These planets are found at orbital distances comparable to the location of 
the boundaries of the habitable zone of their star. 

Aims. In this study, we assess the habitability of planets Gl 581c and Gl 581d (assuming that their actual masses are 
close to their minimum masses) by estimating the locations of the habitable-zone boundaries of the star and discussing 
the uncertainties affecting their determination. An additional purpose of this paper is to provide simplified formulae 
for estimating the edges of the habitable zone. These may be used to evaluate the astrobiological potential of terrestrial 
exoplanets that will hopefully be discovered in the near future. 

Methods. Using results from radiative-convective atmospheric models and constraints from the evolution of Venus and 
Mars, we derive theoretical and empirical habitable distances for stars of F, G, K, and M spectral types. 
Results. Planets Gl 581c and Gl 581d are near to, but outside, what can be considered as the conservative habitable 
zone. Planet 'c' receives 30% more energy from its star than Venus from the Sun, with an increased radiative forcing 
caused by the spectral energy distribution of Gl 581. This planet is thus unlikely to host liquid water, although its 
habitability cannot be positively ruled out by theoretical models due to uncertainties affecting cloud properties and 
cloud cover. Highly reflective clouds covering at least 75% of the day side of the planet could indeed prevent the water 
reservoir from being entirely vaporized. Irradiation conditions of planet 'd' are comparable to those of early Mars, 
which is known to have hosted surface liquid water. Thanks to the greenhouse effect of C02-ice clouds, also invoked to 
explain the early Martian climate, planet 'd' might be a better candidate for the first exoplanet known to be potentially 
habitable. A mixture of several greenhouse gases could also maintain habitable conditions on this planet, although the 
geochemical processes that could stabilize such a super-greenhouse atmosphere are still unknown. 

Key words. Gl 581 - Habitable zone - Darwin - TPF 



1. Introduction 

The M-type star Gl 581 hosts at least 3 planets, which were 
det ected using radial velocity measure ments by Bonfils et 
al. (|2005h (planet 'b') and Udry et al. (|2007ft (planets 'c' 
and 'd'). The properties of this star and its planets are 
given in Table [T] Before this discovery, only two exoplanets 
were known to have a minimum mass below 10 Mg, which 
is usually considered as a boundary between terrestrial and 
giant planets, the latter having a significant fraction of their 
mass in an H2-He envelope. The first one was GJ 876d, a 
very hot planet (P < . 2 da ys) with a minimum mass of 
5.9 M e (jRivera et al.ll2005[ ). The other one is OGLE-05- 
390L b, found to be a ~5.5 Mq cold planet at 2.1 AU from 
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i ts low-mass parent s tar thanks to a microl ensing event 
dBeaulieu et al.l 120061 ; lEhrenreich et all l2006f ). Neither of 
these two planets is considered as habitable, even with very 
loose habitability criteria. In the case of Gl 581, and as 
already mentioned by Udry et al. (2007), the locations of 
planet 'c' and 'd' must be fairly close to the inner and outer 
edges, respectively, of the habitable zone (HZ). In this pa- 
per, we investigate the atmospheric properties that would 
be required to make the habitability of these planets possi- 
ble. 

Because of its equilibrium temperature of ~300 K when 
calculated with an albedo of 0.5, it has been claimed that 
the second planet of this system, Gl 581c, is potentially hab- 
itable (Udry et al. 2007), with climatic conditions possibly 
similar to those prevailing on Earth. After a brief discus- 
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Table 1. Properties of the star Gl 581 and its 3 detected 
planets, from Udry et al. (2007). 



Star 


T efi (K) 


M/M e 


i?/R© 


L/L Q 


Gl 581 


3200 


0.31 


0.38 


0.0135 




Planets 


a (AU) 


M min /M® 


-Rmin/R© 


stellar flux 

S/S *** 


b 


0.041 


15.6 


2.2-2.6 


8.1 


c 


0.073 


5.06 


1.6-2.0 


2.55 


d 


0.253 


8.3 


1.8-2.2 


0.21 



The potential habitability of planets 'c' and 'd', highlighted in 
grey, is discussed in this paper. 

* M m i n = Afsini, where i is the orbital inclination. 

** Radius for a rocky and oc ean planet, respectively l|Sotin et al.l 
llMIValencia et al.|[2007bl ). 
*** So is the solar flux at 1 AU: 1360 W mT 2 . 



sion about the relationship between the equilibrium tem- 
perature and habitability, we summarize in this paper what 
are usually considered as the boundaries of the circumstel- 
lar HZ and the uncertainties on their precise location. In 
Sect. 12.41 we provide parametcrizations to determine such 
limits as a function of the stellar luminosity and effective 
temperature. These can be used to evaluate the potential 
habitability of the terrestrial exoplanets that should soon 
be discovered. We then discuss the specific case of the sys- 
tem around Gl 581. 

2. Habitable planets and the habitable zone 

The HZ is the circumstellar region inside which a terrestrial 
planet can hold permanent liquid water on its surface. A 
terrestrial planet that is found beyond the HZ of its star 
could still harbor life in its subsurface; but being unable to 
use starlight as a source of energy, such endolithic biosphere 
would not be likely to modify its planetary environment in 
an observable way (|Rosingll2005h . In the Solar System, in 
situ searches for biological activity in the subsurface of, for 
instance, Mars or Jupiter's satellite Europa could in prin- 
ciple be carried out. But with exoplanets presumably out 
of reach for in situ exploration, signs of life will have to be 
searched via signatures of photosynthetic processes in the 
spectra of planets found in the HZ of their stars. This is 
the purpose of future space observatories such as Darwin 
(IVolonte et al.ll2000t iKaltenegger fc Fridhmdl I2005D TPF - 



C (jLevine et al.ll2006l ) and TPF-I (jBeichman et al.lll999h . 

For exoplanets, "habitable" thus implies surface habitabil- 
ity. 

A planet found in the HZ is not necessary hab- 
itable. The maintenance of habitable conditions on a 
planet requires various geophysical and geochemical con- 
ditions. Only some of them, those that have a direct 
influence on the atmospheric properti es, are addressed 
in the present paper 



(see for instance Scalo et alj 1200 



Zahnle et all l2ob7jjKasting fc Catlinj 120031 : iLunine) \\M 



Gaidos fc Selsisl 120071 . for a comprehensive view of habit- 
ability). Many factors may prevent (surface) habitability. 
To give several examples: the planet may lack water, the 
rate of large impacts may be too high, the minimum set of 
ingredients necessary for the emergence of life (so far un- 
known) may have not been there, gravity may be too weak 



(as on Mars) to retain a dense atmosphere against escape 
processes and to keep an active geology replenishing the 
atmosphere of CO2, or the planet could have accreted a 
massive H2-HC envelope that would prevent water from be- 
ing liquid by keeping the surface pressure too high. To avoid 
the two last scenarios, the planetary mass should be in the 
approximate range of 0.5-10 M®, although this is more of 
an educated guess than a reliable quantitative estimate. 

Being at the right distance from its star is thus only one 
of the necessary conditions required for a planet to be hab- 
itable. In the current absence of observational constraints, 
we choose to assess the habitable potential of the planets 
with as few hypotheses as possible on their physical and 
chemical nature. We therefore assume that the planet sat- 
isfies only two conditions. Although these two conditions 
are very simple, they may derive from complex geophys- 
ical properties. Future observations will hopefully tell us 
whether such properties are frequent or rare on terrestrial 
exoplanets. These conditions are: 

i) The amount of superficial water must be large enough 
so that the surface can host liquid water for any tem- 
perature between the temperature at the triple point 
of water, 273 K, and the critical temperature of water, 
T c =647 K. This condition implies that the water reser- 
voir produces a surface pressure higher than 220 bars 
when fully vaporized. With an Earth gravity, this cor- 
responds to a 2.2 km layer of water, slightly lower than 
the mean depth of Earth oceans of 2.7 km. For a grav- 
ity twice that of Earth, this pressure corresponds to half 
this depth. Planets with less water may still be habit- 
able, but their HZs may be somewhat narrower than we 
calculate here because liquid water would disappear at 
a lower surface temperature. 

ii) Atmospheric CO2 must accumulate in a planet's atmo- 
sphere whenever the mean surface temperature falls be- 
low 273 K, the freezing point of water. This is a con- 
sequence of the carbonate-silicate cycle, which stabi- 
lizes the long-term surface temperature an d the amount 
of CO 2 in the atmosphere of the Earth (jWalker et al.l 
1 1 9 8 lh . Such an assumption implies that the planet is 
geologically active and continuously outgassing CO2. It 
also implies that carbonates form in the presence of sur- 
face liquid water, which may require continental weath- 
ering. With no atmos pheric CO2, or with a fixed CO2 
level as in Hart (|l979h . the HZ could be ^10 times nar- 
rower than is currently assumed. In the absence of CO2 
(or a greenhouse gas other than H2O), the present Earth 
would be frozen. 

The atmosphere of a habitable planet meeting these condi- 
tions should behave as illustrated in Fig [TJ 



2.1. The equilibrium temperature 

The equilibrium temperature of a planet is given by 



S(l- A ) y 



(1) 



where A is the Bond albedo (which is the fraction of power 
at all wavelengths scattered back out into space - Earth's 
value is 0.29), S is the stellar energy flux, a is the Stefan- 
Boltzmann constant, and / is a redistribution factor. If all 
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distance to the Sun (AU) 

Fig. 1. The CO2 and H2O pressure and the mean surface 
temperature of a habitable planet as a function of the or- 
bital distance. The diagram gives Pco 2 an d Ph 2 (left y- 
axis) and T$ (grey, right y-axis) of an Earth-like planet 
across the inner part of the HZ around the present Sun. 
At orbital distances larger than ~ 0.98 AU, T$ is assumed 
to be fixed at 288 K (its current value on Earth) by the 
carbonate-silicate cycle, which is obviously an idealized pic- 
ture. Beyond 1.3 AU, CO2 condenses in the atmosphere and 
the required level of CO2 depends on the coverage by CO2- 
ice clouds. The pressure of H2O and T s in the inner HZ 
is calculated for a cloud-free atmosphere and 50% cloud 
cover and assuming a reservoir of water that contains (as 
on Earth) more than the equivalent of 220 bars of H2O. 



the incident energy is uniformly distributed on the plan- 
etary sphere, then / = 4. If the energy is uniformly dis- 
tributed over the starlit hemisphere alone, then / = 2. And 
if there is no redistribution, the local equilibrium temper- 
ature can be calculated with / = 1/ cos 9 where 6 is the 
zenith angle. The latter case, for instance, yields good re- 
sults for the surface temperature on the sunlit hemispheres 
of airless bodies with known albedo such as the Moon and 
Mercury. 

It is important to discuss the meaning of T eq and the 
manner in which it can be used to assess habitability. The 
planet Gl 581c has been widely presented as potentially 
habitable because one finds T oq ~ 320 K when calculated 
using the albedo of the Earth. This conclusion is however 
too simplistic for the following two reasons: 

i) For a planet with a dense atmosphere (an inherent prop- 
erty of a habitable planet), T cq does not indicate any 
physical temperatures at the surface or in the atmo- 
sphere. With albedos of 0.75, 0.29, and 0.22, respec- 
tively, and assuming / = 4, Venus, Earth, and Mars 
have equilibrium temperatures of 231 K, 255 K, and 
213 K, while their mean surface temperatures are 737 K, 
288 K and 218 K. The two quantities only match, ap- 
proximately, in the case of Mars, whose tenuous atmo- 
sphere produces a greenhouse warming of only ~ 5 K. 

ii) It can be demonstrated that a necessary (but not suf- 
ficient) condition for habitability is that T eq must be 
lower than about 270 K. If the surface temperature re- 
mains below the critical temperature of water (T c = 
647 K), the thermal emission of a habitable planet can- 
not exceed the runaway greenhouse threshold, ~ 300 W 
m~ 2 (sec Sect. 12.2. 1[) . equivalent to the irradiancc of a 
black-body at 270 K. Therefore, if a planet has an at- 



mosphere and an equilibrium temperature above 270 K, 
two situations may arise. First, T s may remain below T c , 
but there would be no liquid water at the surface and 
no or negligible amounts of water vapor in the atmo- 
sphere. In a second possible situation, the atmosphere 
contains considerable amounts of water vapor, but the 
surface temperature exceeds 1400 K (see Sect. I2.2.1j) . 
This would allow the planet to balance the absorbed 
stellar energy by radiating at visible and radio wave- 
lengths through an atmosphere that is optically thick in 
the infrared (IR). Both cases would render the planet 
uninhabitable. 

For planet Gl 581c to be habitable (i.e., T cq < 270 K), 
its albedo would have to be higher than 0.65. Since planet 
Venus has an albedo of 0.75, this situation may not appear 
unrealistic. However, as we will see in the next sections, 
the Bond albedo is not a quantity given by the planetary 
characteristics alone, but the spectral energy distribution 
of the star also needs to be taken into account. 

2.2. The inner edge of the HZ 

Let us consider a planet with a large water reservoir cover- 
ing its entire surface, but no other greenhouse volatilcs. As 
a first step, we assume that its host star is a Sun-like star 
and that the planet has the same gravity as the Earth. For a 
given orbital distance, a fraction of the water reservoir is in 
the form of vapor. The surface temperature T s imposes the 
surface vapor pressure P w . If P w is high enough, the water 
vapor, in turn, affects T s by blocking the outgoing IR radia- 
tion, by reducing the atmospheric lapse rate and by modify- 
ing the planetary albedo. To account for this coupling, the 
atmospheric structure has to be computed self-consistently 
for a given irradiation. This was done previously by using 
a ID radiative-convective model (Kasting 1988). All the 
orbital distances and stellar fluxes in the following subsec- 
tions are given relative both to the present Sun and to the 
present solar flux at Earth orbit (So = 1360 W m~ 2 ). We 
will see further how these values can be scaled to other stel- 
lar luminosities and effective temperatures. Values given in 
Sect. l2T2TTll2T2~2l and [23^1 were obtained by Kasting ( 1988) 
with a cloud-free radiative convective scheme. Section l2.2.5l 
discusses the likely effects of clouds. 

2.2.1. The runaway greenhouse limit 

For orbital distances smaller than 1 AU (and for the present 
solar luminosity), T s is extremely sensitive to the orbital 
distance, increasing from less than 273 K at 1 AU (in the 
absence of C0 2 ) to about 373 K (P w = 1 bar) at 0.95 AU 
(see Fig. [1]). This sharp increase in T s is mainly caused by 
the increase in the IR opacity and the decrease in the albedo 
caused by absorption of solar near-infrared (NIR) radiation 
by water vapor. For even smaller orbital distances, because 
of the relation between temperature, vapor pressure, and IR 
opacity, the outgoing IR flux becomes nearly independent of 
the surface temperature and tends asymptotically towards 
its upper limit of about 30 W m~ 2 , known as the run- 
away greenhouse threshold (|Abe fc Matsuil Il988t iKastind 
119881: llshiwatari et alj|2002j ). At this point, an increase in 
the irradiation (or a decrease of the orbital distance) does 
not result in an increase in the outgoing IR flux, but leads 
instead to a strong increase in T s and P w . In turn, this 
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produces a slight increase in the albedo and thus in the 
reflected visible/NIR radiation. The increase in the albedo 
for T s above 373 K and P w above 1 bar is the consequence 
of the strong Rayleigh back-scattering occurring in the vis- 
ible (a spectral domain where water vapor does not ab- 
sorb significantly) and to the saturation of the water bands 
absorbing the stellar NIR radiation. This increase in the 
atmospheric albedo, up to about 0.35 (in the absence of 
clouds), protects the water reservoir from complete vapor- 
ization for orbital distances down to 0.84 AU (1.4 So). At 
0.84 AU, T s reaches T c = 647 K, and the water reservoir 
becomes a supercritical fluid envelope. A more limited wa- 
ter reservoir could, of course, be fully vaporized at lower 
irradiation. When this theoretical limit for the irradiation 
is crossed, there is a dramatic increase in T s from > 647 K 
to > 1400 K, a temperature that potentially melts silicates 
on the surface. This behavior is a consequence of the run- 
away greenhouse threshold that limits the mid-IR cooling of 
the planet. At T s > 1400 K, the planet can radiate the ab- 
sorbed stellar energy through the atmosphere at the visible 
and radio wavelengths, at which the water vapor opacity is 
negligible. 

2.2.2. The water loss limit 

In a cloud-free radiative-convective scheme, water vapor 
would become a major atmospheric constituent in an Earth 
analog placed at 0.95 AU from the present Sun. The loss 
of hydrogen to space would no longer be limited by the dif- 
fusion of water vapor from the troposphere to the strato- 
sphere, but by the stellar EUV energy deposited in the 
upper atmosphere, and would be enhanced by ^4 orders 
of magnitude. The hydrogen contained in the whole ter- 
restrial ocean would thus be lost in less than 1 Gyr, which 
would terminate Earth's habitability. For these two reasons, 
0.95 AU could be seen as the inner limit of the present solar 
HZ. 

The water loss limit is difficult to extrapolate to other 
stars and terrestrial exoplanets. Such a limit for present 
Earth corresponds to a surface temperature of ~340 K. 
This is because the background atmospheric pressure is 
1 bar (mainly N2 and O2), and the vapor pressure of H2O 
at 340 K is 0.2 bar, making H2O a 20% constituent of the 
atmosphere. At this H2O mixing ratio, the loss of H to 
space becomes energy-limited. If the background surface 
pressure is higher, then the water-loss limit will be reached 
at a higher temperature. For instance, with a background 
atmospheric pressure of 5 bars, the water loss limit would 
be reached at about 373 K. In addition, the water con- 
tent of terrestrial plan ets is thought to be highly variable 
((Raymond et al.ll2007f ) and could be as high as 50% in mass 
for migrating planets initially formed beyond the snow line 
dKuchnerl 120031: iLeeer et all l2004t I Selsis et all l2007h . For 
such water-rich planets, only a fraction of the water can 
be lost within the lifetime of the planet, and atmospheric 
escape is not a threat to habitability. 

Planets in the Gl 581 system are good ocean planet can- 
didates because the architecture of the system and the high 
mass of the planets are likely to be inherited from type-I 
migration. These planets can thus be initially composed of 
a large fraction of water ice. In addition, the lifetime of the 
water reservoir also depends on the stellar emission in the 
UV, which photodissociates H2O, and in the XUV (0.1- 
100 nm), which produces thcrmospheric heating. The star 



Gl 581 appears quite inactive and should not have strong 
UV or XUV fluxes presently. Thus, gravitational escape is 
not likely to play a major role in the system today. However, 
the high-energy emissions of Gl 581 may have been orders 
of magnitude stronger in the past. The evolution of these 
emissions has not been accurately established for M-type 
stars yet, but they are likely to have had an impact on the 
planet's atmosphere for a more or less extended period of 
time. This point is addressed in more detail in Sect. 13.41 

As already mentioned, the HZ comprises the orbital re- 
gions where terrestrial planets can be probed in search of 
biosignatures. In this context, water loss is an important 
issue, because the leftover oxygen could possibly produce 
a dense Ovrich atmosphere. This may be impossible to 
distinguish from an atmosphere sustained by photosynthe- 
sis, unless additional biomarker gases were also detected. 
Therefore, planets found closer to their parent star than 
the water-loss limit represent questionable astrobiological 
targets. 

2.2.3. The thermophilic limit 

Terrestrial organisms have the extraordinary ability to 
adapt themselves to extreme conditions, including hot envi- 
ronments. For hyperthermophilic prokaryotes, the optimum 
temperature for growth is above 353 K. Pyrolobus fumarii, 
an iron-breathing bacterium, has been found to tolerate 
temperature as high as 394 K, setting the record for the 
highest temperature known to be compatible with life. This 
temperature limit is lower (^333 K) for eukaryotes and pho- 
tosynthetic prokaryotes. The orbital distance at which T s 
cannot be lower than 394 K could thus be seen as an em- 
pirical, anthropocentric, inner edge of the HZ. Interestingly, 
this highest temperature tolerated by life (as we know it) 
is close to the mean surface temperature at which the loss 
of water becomes considerable, so that this limit and the 
water-loss limit are located at approximately the same or- 
bital distance of ^0.94 AU for a cloud-free planet. 

2.2.4. The Venus criterion 

Radar maps of the surface of Venus suggest that liquid wa- 
ter has not been presen t there for at least the last 1 Gyr 
(ISolomon fc Headlll99lD . The Sun was ~8% dimmer at that 
time, accord ing to standard solar evolution models (e.g., 
Baraffe et al. ll998|) . Thus, the solar flux at Venus' orbit was 
then equal to what it is today at a distance of 0.72 AUx 
(1/0. 92) 5 ~0.75 AU. This provides an empirical indica- 
tion of the location of the inner edge of the HZ. We do 
not know if Venus lost its water content after experiencing 
a runaway greenhouse or because its water reservoir was 
much smaller than Earth's. The D/H ratio measured in 
Venusian water vapor traces is 120 ± 40 times higher than 
on Earth (|de Bergh et al.lll991f) . Thus, if the initial D/H 
ratio of Venusian water was the same as on Earth, there 
must have been at least ~120 times more water on early 
Venus than today. This corresponds to an initial inventory 
of about 20 m of precipitable water. Since deuterium is also 
lost to space, albeit at a slower rate than H, and volcanoes 
release poorly deutcratcd hydrogen into Venus' atmosphere, 
this lower limit can significantly underestimate the initial 
water content. This reveals that the cold trap limiting the 
flux of water vapor from the troposphere to the stratosphere 
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did not work on Venus as it worked on Earth. In turn, this 
implies that 1 Gyr ago or longer the mean surface tempera- 
ture was high enough to trigger massive water loss (^340 K 
or higher depending of the background atmospheric pres- 
sure). 

Is it possible to explain the D/H ratio of Venus atmo- 
sphere if we assume that the current C02-rich, hot, and 
mainly dry conditions have always been prevailing at the 
surface of Venus? The answer is yes, if the deuterium en- 
richment has been generated by the loss of the water de- 
livered sporadically by impacts. The minimum amount of 
water lost (~ 10 16 m 3 in volume, 10 19 kg in mass) corre- 
sponds to the contents of 25 000 Halley-sized comets. Venus 
should not have experienced such an accumulated deliv- 
ery since the late heavy bombardment (LHB), which oc- 
curred 3.95-3. 85 Gy r ago. Through numerical simulations, 
Gomes et al. (120051) found that 10 19 kg of cometary mate- 
rial (~ 5 x 10 18 kg of water) impacted on the Moon dur- 
ing the LHB. When scaling this total impact volume to 
Venus (by assuming a factor 1-10 increase relative to the 
Moon) and taking the cometary D/H ratio (twice as large 
as the terrestrial ocean value) into account, we can see that 
enough water could have been brought to Venus during the 
LHB. However, if this water was lost soon after the LHB, 
deuterium has then been escaping during the last 3.85 Gyr 
and only a small fraction should remain today. This implies 
a much higher water delivery that is unlikely to be con- 
sistent with post-accretion impacts. In summary, although 
impact-delivered water is unlikely to explain the present 
D/H enrichment of Venus, we must keep in mind that the 
Venus criterion assumes that the lost water was condensed 
on the surface at some point during Venus' history, which 
has not been proved. 

2.2.5. The effect of clouds on the inner boundary of the HZ 

On a habitable planet close to the inner boundary of the 
HZ, the IR opacity of the atmosphere is fully dominated 
by water vapor, and clouds do not contribute to warm- 
ing the surface, as some types of clouds (high cirrus clouds) 
can do on Earth. But, on the other hand, clouds can signifi- 
cantly increase the planetary albedo and thereby reduce the 
greenhouse warming. In particular, thick clouds forming at 
high altitude, above the level where the incident radiation is 
backscattcred or absorbed, can result in a very high albedo 
and can thus move the habitability limits closer to the star. 
The 1-D radiative convective models that have been used 
to es timate the climatic response to an increase of stellar 
flux (|KastinJl988Lll99lh can only bracket the quantitative 
effects of clouds on the planetary radiation budget. Clouds 
are by nature a 3-D phenomenon that is closely related to 
atmospheric circulation. By adding a cloud layer in a 1-D 
model, it is however possible to investigate the effect of a 
100% cloud cover. We can also estimate the value of T s with 
a 50% cloud coverage (or any other percentage) by assum- 
ing that clouds do not affect the IR outgoing radiation. This 
allows us to calculate the total albedo by combining cloud- 
free and cloudy models and to find the orbital distance at 
which the absorbed energy matches the IR outgoing radi- 
ation. This approximation is acceptable for large enough 
water vapor columns, which are found for T s roughly above 
373 K, and thus for calculations near the inner edge. The 
T s and P w curves for the 50% cloud coverage case in Fig. [1] 



are thus reliable only in the innermost part of the HZ (for 
T s > 373 K). 

A cloud layer located between the 0.1 and 1 bar levels 
(just above the surface in the case of present Earth, but 
at an altitude of ~150 km for a hotter planet with more 
than 100 bar of H2O) produces the maximum increase in 
the albedo. For a planet orbiting our present Sun, these 
clouds can produce an albedo as high as 0.8 and 0.6, re- 
spectively, when covering 100% and 50% of the day side. 
This would move the runaway greenhouse limit to 0.46 AU 
(100%) and 0.68 AU (50%) from the Sun. Interestingly, the 
water-loss limit (T s = 373 K) with 50% cloudiness matches 
the empirical Venus limit (0.72 AU). Such high-altitude and 
thick clouds form mainly by condensation in updrafts, and 
hence are u nlikely to cover a planet's entire hemisphere. 
Joshi (|2003f ) simulated the 3-D climate of a synchronously 
rotating planet (which is expected for planets in circular 
orbits within the HZ of M stars, see Sect. 13. 5j ). with various 
ocean/continent ratios. The author found that the high- 
altitude clouds were mainly located on the dark side of the 
planet and that the cloud cover on the day side was below 
50%. This simulation was calculated for a thin Earth-like 
atmosphere and should not be used to constrain the cloudi- 
ness of a thick H2 0-rich atmosphere close to the inner edge 
of the HZ. It is, however, a good illustration of the simula- 
tions that should be developed in the future to address the 
meteorology of exoplanets. 

2.3. The outer edge of the HZ 

When the carbonate-silicate cycle is at work, the level of 
CO2 should be stabilized at a value sustaining a mean sur- 
face temperature T s somewhere above 273 K. It has to be 
more than a few degrees above 273 K because of the run- 
away ice-albedo feedback that would trigger a global glacial 
event if a relatively large fraction of the oceans is frozen. 
The precise stabilizing temperature should depend on the 
planet's land/ocean distribution, internal heat flow, surface 
carbon inventory, and other factors that are difficult to con- 
strain and that may vary from planet to planet. 

Note that models computing the fluxes between the 
different CO2 reservoirs (atmosphere, ocean, biosphere, 
crust, mantle) and the consequent C O2 atmospheric level 
on an Earth-like p lanet do exist (jFranck et alj l2002t 
ISleep fc Zahnld[2"00lh . Unfortunately, such models typically 
depend on planetary parameters (bulk composition, volatile 
content, radiogenic elements abundances, structure, forma- 
tion history, internal energy, heat flux,...) that are impossi- 
ble to estimate before a planet has actually been discovered 
and deeply characterized. 

Considering the limited understanding of the carbon cy- 
cle on Earth, and the expected diversity of exoplanets, we 
make the very simple assumption that the CO2 level is po- 
tentially maintained at a habitable level. If it is not, then 
the planet is obviously not habitable, which recalls that be- 
ing inside the habitable zone is a necessary but not sufficient 
condition for habitability. 

The relationship between the orbital distance and the 
required level of CO2 is given in Fig.[TJ where the stabiliza- 
tion temperature is arbitrarily fixed to 288 K, the current 
value of T s on Earth. At 1 AU from the present Sun, a 
CO2 abundance of 3 x 10 -4 bar is sufficient to keep T s of 
the Earth at 288 K. Four billion years ago, the solar lumi- 
nosity was about 70% of its present value, and 1000 times 
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more CO2 should have been required to maintain the same 
T s . If the Earth was placed at 1.2 AU from the present 
Sun, the CO2 level would stabilize at approximately this 
same high value. The farthest orbital distance at which the 
Earth could be placed without freezing permanently repre- 
sents the outer boundary of the HZ. This outer edge has 
not been accurately determined yet because of the complex 
effects of clouds resulting from CO2 condensation. 

If we consider a cloud-free CO2 atmosphere (with a 
water pressure fixed by T s ), the outer edge should lie at 
1.67 AU for the present Sun, where a planet would have a 
CO2 pressure of about 8 bar (Kasting 1991). For CO2 pres- 
sures above 6 bar, the cooling caused by the albedo exceeds 
the warming caused by the IR opacity of the CO2 column. 
However, this 1.67 AU limit does not take into account the 
effect of CO2 clouds, which should form a significant cover 
at orbital distances larger than 1.3 AU. 

2.3.1. The C0 2 -cloud limit 

The optical properties of CC>2-ice clouds differ significantly 
from those made of water droplets or H 2 0-ice particles. 
Carbon dioxide clouds are more transparent in the vis- 
ible range, but they efficiently scatter therm al radiation 
around 10 fim. Forget & Pierrehumbert (|l997t ) showed that 
the warming effect caused by the backscattering of the 
IR surface emission exceeds the cooling effect caused by 
the increase of albedo. As a consequence, CO2 conden- 
sation increases the greenhouse warming compared to a 
purely gaseous CO2 atmosphere. Depending on the frac- 
tional cloud cover, the theoretical outer edge of the HZ 
should be found between the 1.67 AU cloud- free limit and a 
100% cloud coyer at 2.4 AU (|Forget fc Pierrehumbertlll997t 
iMischna et al]|2000f) . 

2.3.2. The early Mars criterion 

Numerous geological and geochemical features indicate that 
liquid wate r was present on the surface of Mars as e arly as 
4 Gyr ago (jPollack et al.lll987l ; iBibring et aTll2006f ), when 
the luminosity of the Sun was 28% lower. The solar flux at 
Mars' orbit was then equal to what it is today at an orbital 
distance of 1.5 AUx (1/0. 72) 5 -1.77 AU. Whatever the 
cause of the greenhouse warming on early Mars, this fact 
suggests empirically that the outer edge is located beyond 
this distance. A likely explanation for the early habitabil- 
ity of Mars is the climatic effect of CO2 clouds, perhaps 
combined with additional warming by reduced greenhouse 
gases. Most of the geological features associated with run- 
ning water can indeed be explained by a 3-D model of the 
early Martia n climate, inc luding clouds, with a ~ 1 bar CO2 
atmosphere (|Forgetll2007D . Therefore, the outer boundary 
of the solar HZ must be located somewhere between 1.77 
and 2.4 AU. 

2.3.3. A greenhouse-cocktail limit? 

If a planet's atmosphere contains other greenhouse gases 
in addition to H 2 and CO2, surface conditions could re- 
main habitable at larger orbital distances than the outer 
limit defined by the warming of CO2 clouds. Reduced 
gases, such as methane (CH4) or ammon ia (NH 3 ), could 
for instance increase the greenhouse effect (jSagan fe Chvbal 



Il997f ). Methane could have been an important greenhouse 
gas on early Earth when oxygen was only a trace con- 
stituent of the atmosphere. At that time, the photochemi- 
cal lifetime of atmospheric CH4 was longer and a release of 
this gas at its present rate would have resulted in an atmo- 
spheric level 100-1000 times higher than today. Such a high 
level of CH4 produces a significant greenhouse warming by 
absorbing the thermal IR in the 7.4 /jm band. The efficiency 
of CH4 as a greenhouse gas has, howev er, bee n notably re- 
vised downwards by Haqq-Misra et al. (l2007t ) compared to 
earlier estimates by Pavlov et al. (|2000f ) . The IR absorp- 
tion of CH4 is observable in a global low-resolution mid- 
IR spectrum, as could b e measured by Darwin or TPF-I 
(jKaltenegger et al]|2007t ). In fact, a spectrum of a planet's 
IR emission can potentially reveal any atmospheric species 
working as a major greenhouse gas. 

Since CO2 is not the only greenhouse gas able to main- 
tain habitable conditions at low stellar irradiation, one 
could use atmosphere models to find what the optimum 
composition is to produce a mean surface temperature 
above 273 K at the largest possible orbital distance. This 
could be done by adding species transparent to the stel- 
lar irradiation but whose combined absorption covers the 
whole mid-IR. The combined effect of two or more species 
has to be studied in a self-consistent way. For instance, 
stratospheric warming due to the NIR absorption of CH4 
could induce more cooling than warming by preventing the 
formation of CO2 clouds. Doing the proper modeling, we 
may find that a planet containing a cocktail of, say, CO2, 
CH4, N2O, NH3, and CFCs, may extend surface habitabil- 
ity as far as 3 AU from the Sun (or 0.35 AU from Gl 581, 
see Sect. [3]) The problem with this approach is to justify 
how such a composition could be sustained, taking into ac- 
count the atmospheric sources and sinks, over geological pe- 
riods. Life itself could be a possible answer. A planet could 
be made habitable on purpose by an intelligent civilization 
able to engineer the atmosphere. "Terraforming" of Mars 
has for instance been considered (|McKav et al.lll99~if ). The 
greenhouse cocktail could also be maintained in a "Gaia 
like" homeostasis by some ecosystems that have initially 
originated in the planet subsurface. We will leave this tan- 
talizing problem to further studies and will consider for the 
moment that 2.4 AU is the outermost edge of the HZ for a 
Sun-like star. 

2.4. Limits of the HZ for non solar-type stars 

We have seen that the planetary albedo plays a crucial 
role in the definition of the HZ. When the absorbed stel- 
lar flux per unit surface Sx(l-A)/4 becomes higher than 
300 W m~ 2 , a runaway greenhouse makes the planet un- 
inhabitable. The albedo of the planetary surface is not im- 
portant when determining the inner and outer boundaries 
of the HZ because the albedo of a habitable planet close to 
these edges would be fully determined by its atmospheric 
composition (including clouds) and the spectral distribu- 
tion of the stellar radiation. The Sun emits a large fraction 
of its energy in the visible, a wavelength domain where the 
atmosphere of a habitable planet is highly reflective, be- 
cause of the dependence of Rayleigh backscattering as A~ 4 
and because of the lack of strong H2O absorption bands. 
The emission of stars with a low effective temperature peaks 
in the NIR (around 0.9 /im for Gl 581). In the NIR, the con- 
tribution of Rayleigh backscattering to the albedo becomes 
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Fig. 2. Plots showing the effect of the stellar effective temperature on the albedo. The graphs represent the reflected 
irradiance at the substcllar point of a planet subject to the irradiation from a Sun-like star with T c s = 5800 K (top) and 
from an M-type star with T c s = 3200 K, similar to Gl 581 (bottom). In this illustration, both stars are assumed to radiate 
as black bodies. The graphs give the spectral irradiance of a planet receiving an integrated stellar flux of 1360 W m~ 2 
(left) and 1500 W m -2 (right). The shaded curves with black-body shape give the irradiance of an airless plane t with 
a fixed surface albedo of 0.15. The non-shaded spectra show the irradiance computed with Phoenix JPaillct 20Ql) when 
a cloud-free atmosphere is added. The current Earth atmosphere is shown on the left and a N2-H2O atmosphere with 
T s = 360 K, P s =1.3 bar and P w =0.5 bar on the right. This illustrates the increase in the reflected energy with the stellar 
effective temperature. 



negligible and the strong absorption bands of H2O (plus 
CO2 and possibly CH4 in the outer HZ) cause additional 
absorption of stellar radiation. This effect is illustrated on 
Fig.H 

Because of the relation between the albedo and the ef- 
fective temperature of the star, the limits of the HZ cannot 
be simply scaled to the stellar luminosity. For stellar effec- 
tive te mperatures between 3700 and 7200 K, Kasting et al. 
(|1993t) calculated the albedo of a planet near both edges of 
the HZ that has either a dense H2O or a dense CO2 atmo- 
sphere. Around an M-type star with T e g = 3700 K, for in- 
stance, accounting for this effect represents a 40% difference 
in the stellar flux (15% in orbital distance) for a cloud-free 
H 2 0-rich atmosphere close to the inner edge. The scaling 
factor for a planet with a thin atmosphere, like the modern 
Earth, is smaller - ab out 10% in terms o f stellar flux, or 
5% in semi-major axis (jSegura et ;i 1.1 121)05;. From the vari- 
ous cases studied by Kasting et al. (1993), the limit l- m and 
'out of the solar HZ can be extrapolated to any star with 
luminosity L and an effective temperature between 3700 K 
and 7200 K by using the following relationships: 

'in = ('in© - a^Ti, - b in T?) {^JJ-^j ( 2 ) 



'out — ('out© — Gout?* — b out T^ ^ — j (3) 

with a in = 2.7619 x 10" 5 , b in = 3.8095 x 10" 9 , a out = 
1.3786 x 10~ 4 , 6 out = 1.4286 x 10~ 9 , and T* = T cff - 5700. 
Here, l- m and ' ou t are in AU, and T e g in K. 

As discussed above, the values of '; n Q and 'out© depend 
on the criteria chosen to define the limit of habitability. 
Table[2]givcs the limits of the present solar HZ based on the 



Table 2. Boundaries of the present Solar HZ 





Venus 


clouds 


clouds 


clouds 




crit. 


0% 


50% 


100% 


'in© (AU) 


0.72 


0.84-0.95 


0.68-0.76 


0.46-0.51 




Mars 


clouds 


clouds 


clouds 




crit. 


0% 


50% 


100% 


'out© (AU) 


1.77 


1.67 


1.95 


2.4 



"early Mars" and "recent Venus" criteria, and the radiativc- 
convective models with a cloudiness of 0, 50, and 100%. For 
each cloud coverage, two 'i n © values are given: the runaway 
greenhouse limit, and the T s = 373 K limit. 

Note that the effect of the spectral type on the albedo, 
included in Eqs. (|2|) and {3J as a quadratic function of 
(T e g — 5700), was estimated only for a cloud-free atmo- 
sphere. Since the reflectivity of clouds is less sensitive to 
wavelength, this quadratic term may not be valid to scale 
the boundaries of the HZ for planets covered by clouds. 
The actual boundary locations can be bracketed by using 
the simple luminosity scaling and Eqs. |2]) or ([3]) that in- 
clude the T e ff sensitivity. 

2.4.1. The continuously habitable zone 

The luminosity of a star, and thus the boundaries of its 
HZ, change during its lifetime. Therefore, another criterion 
for planetary habitability could be the time spent at a hab- 
itable distance. For any star of a given spectral type, it 
is possible to define the limits of the region that remains 
within the HZ for a time period longer than a selected 
timcspan. Figure [3] shows the boundaries of the continu- 
ously HZ (CHZ) as a function of the stellar mass, computed 
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Fig. 3. The orbital region that remains continuously hab- 
itable during at least 5 Gyr as a function of the stellar 
mass. The darker area is defined by the empirical "early 
Mars" and "recent Venus" criteria. The light grey region 
gives the theoretical inner (runaway greenhouse) and outer 
limits with 50% cloudiness, with H2O and CO2 clouds, re- 
spectively. The dotted boundaries correspond to the ex- 
treme theoretical limits, found with a 100% cloud cover. 
The dashed line indicates the distance at which a 1 M ffl 
planet on a circular orbit becomes tidally locked in less 
than 1 Gyr. 



with Eqs. @ and © a nd with an evolutio nary model for 
stars of solar metallicity (|Baraffe et a l. 1998). Since Eqs. @ 
and © are not valid for stellar T e g below 3700 K, the effect 
on the albedo was calculated by assuming T e g = 3700 K 
for temperatures below this value. 

3. The cases of Gliese 581c and Gliese 581d 

The star Gl 581 has an inferred effecti ve temperature o f 
3200 K and a luminosity of 0.013 L Q (jUdrv et al.ll2007ft . 
Equations ^ and cannot be applied for T e g < 3700 K 
because no radiative-convective simulation has been per- 
formed for such low stellar effective temperatures. However, 
the albedo calculated for T cS = 3700 K and 3200 K should 
be similar. Figure 0] shows the location of the boundaries of 
the HZ around Gl 581 for various criteria, along with the 
three known planets. The limits are given with or without 
the albedo correction (except for the cloud-free simulation, 
which was specifically run for T c s — 3700 K). 

3.1. The planet Gl 581c 

Figure [4] shows that planet Gl 581c would be habitable 
only if clouds with the highest reflectivity covered most of 
the daytime hemisphere. A 50% cloud cover is not enough 
to prevent a runaway greenhouse effect on Gl 581c, which 
receives 30% more energy than Venus today. This problem 
is exacerbated by the fact that Venus has a much higher 
albedo than the expected value for a habitable planet at 
the orbital distance of Gl 581c. The composition of the 
atmosphere of Gl 581c depends on the mass of the initial 
water reservoir on the planet, and on the efficiency of the 
gravitational escape of H. Two possible scenarios, inspired 



by the fate of Venus and by simulations done for an Earth- 
like planets, can be suggested: 

i) The reservoir of water is large so only a fraction of it 
will be lost to space. Numerical simulations suggest that 
planets more massive than about 5 M m , such as thos e 
found around Gl 581 or HD 69830 (jAlibert et al.ll2006f ). 
started their formation in the cold outer protoplanetary 
disk, accreting some icy planetesimals, and migrated 
close to the star. In this situation, water can be orders of 
magnitude more abundant than on Earth. If this is the 
case for Gl 581c, its rocky surface is probably covered by 
a thick layer of H2O. Depending on the amount, water 
can form a mantle of hot and high-pressure ice under- 
neath a fluid envelope of supercritical H2O (Selsis et al. 
2007). 

ii) The reservoir of H2O has been lost earlier, when the 
star was still a strong X-ray and EUV emitter. A CO2- 
rich, Venus-like atmosphere could remain, with surface 
temperatures still in the range 700-1000 K. A significant 
level of leftover O2 may be present. 

But the planet could also have a bulk composition that is 
so different from Earth that a broader diversity of situa- 
tions should be explo red. The planet could, for instanc e, be 
a volatile-rich planet (lKuchnerll2003tlLeger et al.ll2004D or a 
carbo n- or ceramic-planet ( jKuchner fc Sc agcr 20051; iGaidod 
|2000D , with an atmospheric composition very different from 
our assumptions derived from the study of Venus, Earth, 
and Mars. However, water is considered an indispensable 
prerequisite for life and any additional atmospheric con- 
stituent would likely enhance the opacity of the atmosphere 
and the amount of greenhouse warming. The inner edge de- 
termined for a pure H2O atmosphere thus seems to be the 
closest possible limit for habitability. 

3.1.1. The effect of gravity 

The gravity at the surface of Gl 581c, calculated for its min- 
imum mass and for the range of radii predicted by Sotin et 
al. (2007) or Valencia et al. (2007), is between 1.3 and 2 g. 
Therefore, a given value of the vapor pressure P W (T S ) would 
correspond to a vapor column 1.3-2.0 times smaller than 
the equivalent at Earth gravity. As the water column de- 
termines the IR opacity of the atmosphere, a higher gravity 
is expected to shift the inner edge of the HZ closer to the 
star. However, IR opacity is not the only parameter to be 
affected by the gravity. For a given value of T s , gravity also 
influences the lapse rate and the albedo produced by the 
water column. When treated self-consistcntly, these effects 
tend to compensate for each other, thus weakening the over- 
all influence of gravity. Radiative-convective simulations of 
a planet with a s urface gravity of 2.5 g were carried out 
by Kasting et al. (|l993f ). who found an inner edge only 3% 
closer. For the expected gravity of Gl 581c, its effect can be 
safely neglected. 

3.1.2. Stability of a quasi-snowball state 

For a given stellar irradiation, there is not necessarily a 
unique solution for T s . On Earth, for instance, climate is 
known to be bistable. It has been warm during most of 
its history, with a mean surface temperature well above 
273 K (which includes glacial and interglacial periods), ex- 
cept for a few snowball Earth events characterized by an 
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Fig. 4. Diagrams depicting the HZ around the Sun and Gliese 581. The grey areas indicate the theoretical inner edge 
for different fractional cloud covers. The width of each inner edge is defined by the runaway greenhouse and water loss 
limits. The thick lines give the inner and outer boundaries of the "empirical" HZ, based on the non-habitability of Venus 
for the last 1 Gyr and the apparent habitability of Mars 4 Gyr ago. The dashed line gives the outermost theoretical limit 
of the HZ, found with a 100% CO2 cloud cover. The upper diagram shows the limits computed for the Sun's properties. 
The lower diagram shows the limits computed for a 3700 K M-type star and scaled to the stellar luminosity of Gl 581 
(for which T e s = 3200 K). The limits marked with the symbol * are simply scaled to the luminosity and are thus slightly 
closer to the star - they are shown here because the wavelength dependency of the albedo has been explicitly computed 
only for thick and cloud- free H2O- and C02-rich atmospheres. The albedo should be less sensitive to the stellar effective 
temperature in the presence of clouds. For the eccentric solution found by Udry et al. (2007), the orbits are confined 
within the grey bars about the marked planet locations. 



ice cover down to the equator. The most recent snowball 
events occurred at the end of Ncoprotcrozoic era between 
730 and 610 Myr ago. At the beginning of these events, 
the runaway ice-albedo feedback makes the global mean 
temperature drop to ^220 K for a few tens of thousands 
of years. This temperature drop is followed by a period 
of a fe w million years wi th a mean temperature around 
265 K (|Schrag et al.ll2002T ). Thanks to the volcanic release 
of CO2 and the inefficiency of carbonate formation on a 
frozen Earth, our planet did not remain trapped in such 
state. This constitutes strong evidence of the long-term sta- 
bilization of the climate through the carbonate-silicate cy- 
cle. Snowball events illustrate the fact that the Earth would 
not be close enough to the Sun to maintain T s above 273 K 
if water vapor was the unique greenhouse gas in the atmo- 
sphere. Without enough atmospheric CO2, the surface of 
our planet would be frozen, and the very high albedo of the 
frozen oceans would permanently keep T s as low as 220 K. 

It can be considered whether such bistability could be 
possible on Gl 581c. Radiative-convective models give one 
solution for the atmosphere structure and surface tempera- 
ture, but would a cold solution be possible? Let us assume 
that ice covers the entire planetary surface, resulting in an 
albedo of ~0.8, and that the atmosphere has a negligible 
radiative effect, either because it is very tenuous or because 
it is composed of IR-transparent gases like N 2 . In this case, 
T eq would be close to 235 K, which is consistent with this 
"snowball" hypothesis. However, in the absence of a green- 
house effect, and in the case of a slowly rotating planet 
(see Sect. I3.5p . the surface temperature at the substellar 
point would be 330 K. This temperature corresponds to a 
vapor pressure of 0.2 bar, which is inconsistent with the as- 



sumed albedo and atmospheric transparency. If we assume 
an extreme albedo of 0.95, as on Enccladus, the icy satellite 
of Saturn, the substellar temperature would be well below 
273 K, and the icy state would be stable. However, in the 
absence of liquid water on the surface, the planet would not 
be habitable. An interesting situation might be found for 
a very narrow range of albedo (^4 = 0.89 — 0.90) for which 
most of the planetary surface is frozen, except for an area 
at the substellar point, where the local equilibrium tem- 
perature and the vapor pressure slightly exceed 273 K and 
6.1 mbar, respectively. This would allow for the presence of 
liquid water. 

In this latter case, though, the IR opacity and effect on 
the albedo of a 6.1 mbar vapor column is already impor- 
tant (in fact the mean water content of Earth's atmosphere 
is very close to this value) and is not consistent with the 
assumed surface temperature. This means that tempera- 
ture and vapor pressure are expected to diverge rapidly 
from this assumed starting point. It would be interesting 
to study this case with a time-dependent model includ- 
ing the vaporization of water at the substellar point and 
its condensation in the cold regions, to check whether a 
substellar hot spot can exist without triggering a runaway 
greenhouse. Such a model should also include the rotation 
of the planet and the variation in the orbital distance due to 
the eccentricity. Indeed, a partially habitable steady-state 
implies that the substellar point moves on the planetary 
surface (otherwise all the water inventory ends up in the 
frozen regions), which requires significant eccentricity (see 
Sect. 13 . 5]) . However, this case is ad hoc as there is no partic- 
ular reason why the planetary albedo would have this exact 
value and, if the planet was mainly frozen, any volcanic re- 
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lease of CO2 (expected on a planet as massive as Gl 581c) 
would trigger a runaway greenhouse. 

3.2. The planet Gl 581d 

Planet Gl 581d receives about half the energy flux that 
Mars gets currently from the Sun. However, because CO2- 
rich atmospheres absorb more energy from an M-type star 
than from a G-type star, the orbital distance of Gl 581c is 
no more than 4% beyond the empirical "early Mars" limit 
(when scaled to Gl 581). Moreover, a partial cover of CO2- 
ice clouds can theoretically sustain habitable conditions for 
even lower stellar fluxes. The third planet of this system is 
thus potentially habitable, according to our present knowl- 
edge. Its high mass makes the maintenance of a thick at- 
mosphere possible over billions of years. Of course, the high 
(~ 8 M®) minimum mass also means that the actual mass 
could be > 10 M®, and the planet may thus be a gas or ice 
giant, rather than a rocky planet. If so, this planet is only 
one of the many other gas giant exoplanets already known 
to be within the HZs of their parent stars. 

The warming effect of C02-ice clouds has only been 
studied for Sun-like irradiation. Because of their size and 
their optical constants, typical C02-ice particles are rela- 
tively transparent to visible radiation, but scatter efficiently 
at mid-IR wavelengths. For a stellar flux shifted towards the 
NIR, the albedo of these clouds could be significantly higher 
than for solar irradiation. This effect is not included in 
Eq. Detailed modeling is thus required, without which 
it is not possible to determine the precise location of the 
edge of the HZ as defined from C02-ice clouds, and thus the 
potential habitability of the large terrestrial planet Gl 581d. 

As discussed in Sect. 12.3.31 some atmospheric compo- 
sitions with the right greenhouse cocktail might provide 
stronger greenhouse warming than CO2 alone and sustain 
liquid water at the surface of Gl 581d. In addition, some 
greenhouse gases, such as CH4, condense at much lower 
temperatures than CO2 (the triple point of CH4 is at 90 K 
and 0.1 bar) and would not be trapped as ice on the night 
side in case of synchronous rotation (see Sect. 13.5)) . It is not 
obvious that CH4 alone could maintain habitable conditions 
on Gl 581d. Indeed, the efficiency of the surface warming 
is expected to decrease above a certain CH4 level. This is 
due to the stratospheric heating by the direct absorption 
of the visible and near-IR stellar flux. For Sun-like irradi- 
ation and partial pressures of CH4 above about 5 mbars, 
this effect competes with the surface greenhouse warming. 
The temperature structure of a CH4-rich atmosphere under 
the irradiation of the star Gl 581 should thus be specifically 
calculated to answer this question. A remaining problem is 
to identify the geochemical mechanisms able to stabilize a 
CH4-rich or a "super-greenhouse" atmosphere over several 
Gyr. 



3.3. An excess of volatiles? 

The volatiles that constitute the atmosphere of terrestrial 
planets in the Solar System have been mainly accreted 
as solids (rocks, with a possible minor contribution from 
ices). This is inferre d from the c omposition in noble gases 
and their isotopes (|Pepinl 1991). If a significant mass of 
H2 and He was accreted directly from the protoplanetary 
nebula, most of it was lost to space and the traces that 



remain represent a negligible fraction of the present at- 
mospheres. The last phases of the accretion of Earth are 
likely to have occurred after the dissipation of the protoso- 
lar nebula. Indeed, disks are fo und to typically live less than 
5-10 Myr (jMever et alj|2007t) . while it shoul d have taken 
more than about 30 Myr to form the Earth (|Martin et al.l 
120061 ). However, due to the diversity of disk properties (in 
lifetime, mass, density) and the role of migration (which ac- 
celerates the accretion), it is possible that the accretion of 
nebular gas is much more efficient on some terrestrial proto- 
exoplane ts, esp ecially if they are more massive than Earth. 
Rafikov (2006) found that planets more massive than 6 M© 
could have accreted a significant fraction of their mass as 
H2-He. For instance, a 8 M© planet could be made of 7 M© 
of rocks and 1 Mq of H2-He. The dissipation of the disk 
could have frustrated the evolution of such planet into a gas 
giant. It is thus possible that some planets more massive 
than ~6 M ffi (which is potentially the case of Gl 581c and 
Gl 581 d) could not host liquid water because of the pres- 
sure and temperature imposed by a massive envelope of gas. 
This situation would be similar to that of HP 69830d. This 
18 M e planet found in the HZ of its star (jLovis et al.ll2006h 
is thought to have accreted a massive H2-HC envelope, un- 
derneath which wate r can only exist as a supercritical fluid 
or high-pressure ice (jAlibert et al.ll2006D . 

As mentioned before, the Gl 581 planets may very well 
have started their formation beyond the snow line. In this 
case, and even without invoking the accretion of hydrogen- 
rich gas from the protoplanetary nebula, their volatile con- 
tent could be orders of magnitude higher than on Earth. 
If only 10% of the accreted solids are made of cometary 
ice, a 6 M e planet would contain as much as 0.55 M ffi of 
water and 0.06 M© of other volatiles, mainly CO2, CO, 
CH4, and NH3. For comparison the mass of the terrestrial 
oceans is 2 x 10" 4 M ffi and the mass of Earth's atmosphere 
is 9 x 10~ 7 M®. 

The distribution of these volatiles between the interior 
and the atmosphere (including here surface ices) at the end 
of the accretion depends on the planet's thermal history. 
Also, the volatiles initially present in the planetesimals may 
not be retaine d during the violent accrctional collisions 
(jLissa ucr 2003). However, the planet can possibly have a 
massive and enriched gaseous envelope preventing the liq- 
uid phase of water from existing on the surface. 



3.4. Atmospheric erosion and the habitability of planets 
around M-type stars 

Models presented in this paper assume that planets can 
retain a dense atmosphere. However, planets in the HZ of 
active stars can be exposed to high levels of X-ray and EUV 
radiation (XUV) and strong particle fluxes from the quies- 
cent stellar wind or coronal mass ejections (CMEs). Such 
high emissions are the result of the stellar magnetic activity 
and can induce important thermal and non-thermal atmo- 
spheric losse s to space, potentially able to strip the whole 
atmosph e re (|Griefimeier et al.l 120051 : Lammer et al.l 120061: 

^^h| I 



LammeJ 
20071: IScalo et al 



2007T iKhodachenko et al.l 120071 : I Lammer et al.l 



2007;). Within the HZ of the Sun and 



solar-type stars, conditions threatening the survival of the 
atmosphere and the ha bitability are limited to the first few 
hundred million years (|Kulikov et al.ll2007f) . and we know 
that the Earth and Venus atmospheres survived this early 
active phase (more damage may have been caused to Mars) . 
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But extreme irradiation conditions could last several Gyr in 
the case of M-type stars, as discussed below. In the context 
of Gl 581, it is thus worth addressing the question of stellar 
activity and the impact on the habitability of its planets. 

For an initial estimate of the evolution of XUV irra- 
diances we have used a proxy indicator, which is the ra- 
tio of the X-ray luminosity to the bolometric luminosity 
(log [Lx/iboi])- This ratio is highest for the most active 
stars (i.e. fastest rotation period) and decreases monotoni- 
cally with decreasing l evel o f chromospheric activit y (e.g., 
Stelzer & Neuhauser, 120011 : Pizzolato et al. l2003h . From 
the analysis of open-cluster stars it is now well estab- 
lished that all single late-type stars (G-K-M) spin down as 
they age. Their activity decreases with time, and so does 
the ratio log(Lx/^boi)- It is also a well-known effect that 
log(Lx/^boi) docs not increase up to values arbitrarily close 
to unity for very active stars. Instead, a "saturation" phe- 
nomenon occurs and no active star seems to go higher (ex- 
cept for the stro ngest flares) than log( Lx/L boi) ~ —3 (e.g., 
Vilhu & Walter. KL987I ; Stauffer et al.. ll994h . Qualitatively, 
the evolution of log(Lx/^boi) for a late-type star has a 
flat plateau from its arrival on the main sequence up to a 
certain age (end of saturation phase) and then decreases 
monotonically as a power law function of age. 

We have compiled a sample of K- and M-type stars with 
ages determined in a similar way to the solar analogs in the 
"Sun in Time" sample (Ribas et al. 2005), i.e., memb ership 
in clusters and moving groups (jMontes et al.ll200lh . rota- 
tion period, membership in wide binaries, and isochrones. 
For these, log(Lx/£boi) values have been obtained from 
the thorough list provided by Pizzolato et al. (2003) and 
complemented in a few cases with values estimated directl y 
from ROSAT measurements following Schmitt et al. ljl995f ). 
The evolution of log(Lx/^boi) with age for stars of differ- 
ent spectral types is illustrated in Fig. [5] Note that today's 
Sun has a value of log(Lx/Lboi) ~ —6.1. The solid lines 
are semi-empirical estimates that result from the fact that 
the evolution of Lx in the power law regime is roughly in- 
dependent of the spectral type and Lboi comes from the 
use of stellar models for different masses. For the K- and 
M-type stars these are preliminary estimates and the un- 
certainty of each point can be of a few tenths of a dcx. In 
the case of G-type stars, the values plotted are much more 
reliable and come from direct measurements for the "Sun 
in Time" targets. More details on this analysis will be given 
in a forthcoming publication (Ribas et al. 2008, in prep.). 

It is reasonable to assume that the overall XUV flux 
received by a planet in the HZ scales with the ratio 
log(Lxuv / ^boi) ■ Figure [5] shows that, while solar-type stars 
stay at saturated emission levels until ages of ~100 Myr 
and then their XUV luminosities rapidly decrease follow- 
ing a power-law relationship as a function of age, M-type 
stars have saturated emissions (i.e., highest activity) up to 
ayes of a few Gyr. Thus, planets in the HZ of M-type stars 
may receive XUV fluxes that are 10-100 times higher than 
in the HZ of solar-type stars of the same age. 

To estimate the possible activity level of Gl 581 we 
carried out a sear ch in the ROSAT Faint Source Catalog 
(|Voges et al.ll2000f ) that produced no result. Thus, Gl 581's 
X-ray flux is below the ROSAT detection threshold. Such a 
threshold can be calculated following Schmitt et al. (1995) 
and we found that Lx must be lower than 10 27 erg s _1 at 
the distance of Gl 581, implying that log(Lx/Lboi) < —4.7. 
From Fig. O this leads to a lower limit of the age that, 




-0.5 0.0 

log(t [Gyr]) 

Fig. 5. Evolution of the ratio between the X-ray and bolo- 
metric luminosities as a function of age for stars of different 
masses. The solid lines represent semi-empirical laws, while 
symbols give observed values for G (+), K (*) and M (o) 
stars. The dashed line gives the upper limit for the value 
of log(Lx/Lboi) in the case of Gl 581, for which there is no 
ROSAT detection. 



considering the associated uncertainties, could be around 
7 Gyr. However, we caution that this is still a prelim- 
inary estimate and that Gl 581 should be studied fur- 
ther before giving a conclusive figure. With regards to 
an upper limit to the age, the space motions of Gl 581 
(U, V, W) = (-22.7, -29.3, -8.0) km s _1 are characteristic 
of a galactic disk star, and thus its age should not be above 
10-11 Gyr. Also, the star is known to have a only midly 
subsolar metalli city ([Fe/H m -0.3; Bonfils et al. 120051 : 
Bean et al. l2006l ). which is also consistent with being a disk 
member. 

The current XUV emissions of Gl 581 may not represent 
a significant threat to the stability of planetary atmospheres 
inside the HZ, but these planets have most likely suffered 
strong erosion phenomena during earlier stages. Dominant 
loss mechanisms are thermal escape, induced by the heat- 
ing of the upper atmospheric layers by XUV irradiation, 
and non-thermal escape caused by the interaction of these 
upper layers with the stellar particle wind. The efficiency 
of these erosion processes depends on several factors. First, 
the composition of the atmosphere controls the tempera- 
ture of the exosphere and thus the thermal escape rates, 
which ca n be extremely high if the blow-off temperature is 
reached (jLammer et al.ll2"006|) . A high mixing ratio of CO2, 
for instance, allows the exosphere to cool efficiently by IR 
emission (that is why the exosphere of Venus is about three 
times cooler than Earth's). The thermal structure of the 
upper atmosphere also determines the extension of the ex- 
osphere and indirectly the efficiency of non-thermal losses 
(like sputtering, ion-pick up, dissociative recombination). 
Non-thermal losses also depend on poorly constrained pa- 
rameters like the intrinsic magnetic moment of the planet 
and on the stellar wind properties (velocity and density). 

Observational constraints on stellar winds exist 
( Wood et al.112005 1 ). but the ir interpretation is still debated 
( Holzwarth fc JardineH2007[ ). It does however seem well es- 
tablished that the particle winds of active M-type stars 
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could at least be 10 times stronger than that of today's Sun, 
in the HZ. During the early and active stages of an M-type 
star like Gl 581, the rate of atmospheric losses due to the 
association of XUV heating and particle wind-induced ion 
pick-up could theoretically reach tens of bars of heavy ele- 
ments (O, C, N) per Gyr for a weakly-magn etized planet in 
the H Z, even with a COvrich atmosphere (jLammer et al.l 
120071) . Tidal interactions rapidly slow down the rotation of 
planets in the HZ of M-type stars (see Sect. I3.5[> . which 
may result in weak intrinsic magnetic moments, unable to 
provide sufficient protection against non-thermal erosion. 
On the other hand, Gl 581c and Gl 581d are significantly 
more massive than the Earth, which may result in an en- 
hanced dynamo and better shielding. Their surface gravities 
are also higher by factors 1.3 and 2.5, respectively, which 
limits the extension of the exosphere and should lower the 
erosion rate. 

Eventually, the threat posed by the atmospheric ero- 
sion depends on the inventory of volatiles, as has already 
been discussed in the case of the water loss in Sect. 12.2.21 
If the loss rate exceeds the outgassing rate, or if the at- 
mospheric losses occurring during the planet's early stages 
exhaust species essential for the maintenance of life, the 
planet would obviously not r emain habitable in the long 
term. Gaidos & Selsis (|2007f ) point out that a replica of 
early Earth (assumed to have a C02-rich atmosphere) lo- 
cated in the HZ of a M-type star could lose its entire at- 
mospheric nitrogen inventory during its first Gyr via ther- 
mal escape from XUV heating alone. A large reservoir of 
volatiles and/or a continuous outgassing of the essential 
atmospheric species is thus required to ensure durable hab- 
itable conditions around M-type stars. 

Models providing quantitative estimates for atmo- 
spheric escape in the extreme conditions encountered in the 
HZ of M stars are still under development. Available prelim- 
inary results cited in this section are extremely promising 
and stress the importance of atmospheric retention for hab- 
itability around low-mass stars. However, current estimates 
of thermal and non-thermal losses in this context assume 
a static upper atmosphere, while high atmospheric escape 
rates would require a hydrodynamic treatment. This field of 
research has recently become very active due to the impor- 
tant questions raised by the habitability of planets around 
M-type stars. Although we have to end this section with 
a question mark, the reader should expect forthcoming re- 
sults obtained with thermal hydrodynamic models coupled 
with non-thermal erosion that will shed more light on the 
fate of planetary atmospheres around active stars. 

3.5. The effect of rotation and eccentricity 

Udry et al. (2007) present two fits of the radial velocity 
curve, one obtained with fixed circular orbits, and another 
where the eccentricity is a free parameter. The latter re- 
sulted in the best agreement and yielded e = 0.16 and 
e = 0.2 for Gl 581c and Gl 581d, respectively. 

If the orbits are circular, then these planets could have 
evolved into a synchronized state, in which the rotation 
period (P r ) equals the orbital period (P). The synchroniza- 
tion time of these planets is indeed very short. We obtain 
25 000 yr and 10 7 yr for planets Gl 581c and Gl 581d, re- 
spectively, when adopting a value of Q = 100 for the dissi- 
pation factor, typical of terrestrial planets. In this case, the 
planets have permanent day and night hemispheres. It is in- 



teresting to note that we find a much longer circularization 
time for these two planets (10 10 yr), which is consistent 
with eccentric orbits. If the planets do have a significant 
eccentricity, then the following cases are possible: 

i) Because terrestrial planets are expected to have a non 
purely axisymmetric shape, tidal spin down may have 
made it possible for the rotation to be captured in a 
spin-orbit resonance with P = (n/2)P r . Mercury, for 
instance, is trapped in a 3:2 spin-orbit resonance as a 
consequence of its oblong shape. For resonances with 
n 7^ 2, the star rises and sets at any point of the planet 
and there is no permanent dark area. In the synchronous 
case (n = 2) , there is a region of the planet that never re- 
ceives direct starlight. For an eccentric planet with zero 
obliquity, the size of the corresponding surface is gov- 
erned by the amplitude of the optical librations about 
the substcllar point when neglecting the smaller contri- 
bution of forced and free librations. Considering up to 
third-order terms in eccentricity, we estimate that the 
amplitude of the libration is 2 e + 25 e 3 /16, and thus the 
fractional area r of the planet that can receive direct 
starlight, is then given by 

1 2e + 25e 3 /16 
r = x + — • (4) 

For e = 0.2, the optical libration amounts to ~ 23°, 
and r would take a value of 63%. Any non-zero obliq- 
uity adds latitudinal librations that may significantly 
increase these values. For high obliquities, it is possi- 
ble for the entire surface of the planet to receive direct 
starlight some time during an orbital period. However, 
high obliquities are not consistent with tidally-driven 
evolution models of the spin axis. 
ii) The planet has not been captured in a spin-orbit res- 
onance. In this case, no region of the planet is perma- 
nently dark. 

Let us now consider the synchronized case of P r = P, 
with low obliquity. Whether the orbit is eccentric or not, 
there is a region on the planet that is not directly heated by 
the stellar radiation. Such a region can be an entire hemi- 
sphere if the eccentricity is negligible. There is thus a risk of 
the atmospheric constituent condensing irreversibly on an 
ice sheet covering the dark region. Joshi et al. (1998) and 
Joshi (2003) have shown that the atmosphere has to be suf- 
ficiently dense and opaque in the IR to transport enough 
heat to the dark area and to prevent such an atmospheric 
collapse. Near the edges of the HZ, a habitable planet ful- 
fills these conditions because it is very rich in either H2O 
or CO2. Therefore, if one concludes that planets Gl 581c 
and Gl 581d are habitable thanks to the effect of clouds, 
then synchronized rotation would not appear to represent 
an impossible obstacle to their habitability (although this 
would give rise to a climate that is certainly not "Earth 
like", with a super-rotating atmosphere). In addition, the 
strong irradiation of planet Gl 581c is expected to produce 
strong thermal tides in the atmosphere (similar to those on 
Venus) that could preven t the synchronization of the planet 
(|Correia fc Laskarfl2003h . 

For planet Gl 581d, whose orbit is in the outer area of 
the HZ, synchronization can pose a more significant threat. 
The habitability of the surface of the Earth has remained 
constant during most of its history, except for a few snow- 
ball events possibly generated by climatic instabilities. The 
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Earth was able to recover from these frozen episodes thanks 
to the continuous release of volcanic CC>2- However, a syn- 
chronized planet may not be able to recover from a snowball 
event as effectively as the Earth did because the conden- 
sation of CO2 on the dark side would occur much faster 
than the release of volcanic CC>2- Only a rapid and dra- 
matic event, such as a large impact on the ice sheet, could 
initiate a new habitable period. The habitable state of a 
synchronized planet could thus be particularly fragile when 
it relies on the presence of atmospheric CO2, which is the 
case for most of the HZ except for the very narrow inner 
part where H2O is "self sufficient". 

Long-term habitability on a synchronized planet could 
still be maintained if the greenhouse warming was provided 
by atmospheric compounds that remain gaseous at suffi- 
ciently low temperature, such as CH4. In this situation, 
one may wonder if an inefficient transport of the incident 
energy from the starlit to the dark hemisphere could help 
maintain habitable conditions on the starlit hemisphere, or 
in a smaller region around the substellar point, while the 
dark hemisphere would remain too cold to host liquid wa- 
ter. In this case, restricting the habitability to a fraction of 
the planetary surface would compensate for the low stellar 
flux. Although this case represents a tantalizing configura- 
tion, the cold trap for water would still operate, and this 
would be irreversibly transported from the habitable to the 
frozen region. 

In the case of the planets having an eccentricity as high 
as 0.16-0.2, the study of the climate at steady state is of 
course less relevant than for a circular orbit. As shown 
by Fig. g both Gl 581c and Gl 581d are likely to make 
incursions inside the HZ (or excursions outside the HZ) 
depending on the assumed boundaries. This situation is 
not expected to change the discussion significantly regard- 
ing planet Gl 581c, since the timcscales associated with 
the vaporization of the water reservoir are significantly 
longer than the 12-d orbital period. In addition, Williams 
& Pollard (|2002h show that habitability depends primar- 
ily on the average stellar flux received over an entire orbit, 
< S >, even at high eccentricity, and < S > depends on 
the eccentricity through the relationship: 

where S is the stellar flux at the semi-major axis of the 
orbit. For e = 0.16 and e = 0.20 the average stellar flux is 
greater by 1.3% and 2.1%, respectively, in comparison with 
a circular orbit. With Gl 581d located near the outer edge 
of the HZ, such slightly higher flux may be precious for the 
maintenance of habitability. Excursions out of the HZ may, 
however, increase the efficiency of irreversible condensation 
on the dark side, although these excursions could last only 
a fraction of the 83-d period. 

3.6. Internal heat flux and tidal dissipation 

Geothermal heat has a major impact on climate during 
the early phases of planetary evolution. Both during and 
shortly after accretion, it can trigger a runaway green- 
house and a surface magma ocean provided that the cool- 
ing to space is limited by a dense at mosphere, which woul d 
be blanketing the thermal emission (|Abe fc Matsuilll988h . 
This is due to the large amount of energy released in a short 



period of time and to an efficient convective transport of in- 
ternal energy to the surface by a yet liquid mantle. The pe- 
riod during which Earth's surface was significantly heated 
"from underneath" was restricted to the first tens of Myr. 
On the present Earth, the internal heat flux, dominated by 
the decay of radiogenic species and the release of the initial 
accretion energy, is about 5000 times lower than the solar 
energy absorbed by the planet. Thus, the direct influence of 
the internal heat on the climate has been negligible for the 
past 4.50 Gyrs or so. Similarly, the Gl 581 system is most 
likely too old (see Sect. I3.4[) to invoke geothermal forcing 
of the climate. 

The difference in mass and radius between the Earth 
and the Gl 581 planets should not affect this reasoning sig- 
nificantly. The mean flux of accretion energy through the 
surface scales with M 2 /R 3 and remains within 10 times 
that of the Earth. Assuming the same abundance of radio- 
genic species for the Gl 581 planets, the associated heat 
flux should scale with M/R 2 and remain Earth-like within 
a factor of ~2. Only by assuming orders of magnitude more 
radiogenic species per unit mass could the geothermal flux 
be significant for the global energy balance on the surface. 

Interestingly, during the few Myr after the Moon- 
forming impact, the tidal interaction with the Moon might 
have been a majo r heating source for Earth's climate 
(jZahnle et al.ll2007t ). Tidal dissipation was favored by the 
small Earth-Moon distance and most of the dissipated en- 
ergy was taken from the kinetic rotation energy. Because 
the rotation of Gl 581 planets is presumed to be tidally 
evolved, tidal dissipation may still occur but through the 
damping of their high eccentricity and may increase the in- 
ternal heat flux to a level that could affect the climate. It 
is thus worth exploring this point further. 

Assuming synchronously rota ting a nd eccentric planets, 
and using formulae from Peale (|1999T ) with a tidal dissi- 
pation factor Q = 100, we find that 1.4 x 10 16 W and 
3.1 x 10 12 W are dissipated respectively in planets Gl 581c 
and Gl 581d (these values are obtained with the largest radii 
given in Table [T] to analyze the maximum effect). This dis- 
sipated power corresponds to a maximum surface flux of 
10 W m~ 2 for planet Gl 581c and 0.0015 W m~ 2 for planet 
Gl 581d. 

On Gl 581c, tidal dissipation can contribute slightly to 
the energy balance, with the heat flux being potentially as 
high as 2% of the absorbed stellar energy (about 1% of 
the stellar flux at the planet's distance, averaged over the 
planetary surface). Besides this direct warming, this heat 
flux can increase the volcanic activity and the release of 
CO2. At the orbital distance of Gl 581c, additional CO2 can 
only further destabilize the climate against water loss or a 
runaway greenhouse. Therefore, tidal heating caused by a 
possibly non-zero eccentricity speaks against any habitable 
conditions on this planet. 

On planet Gl 581d, tidal dissipation is negligible 
(0.003% of the absorbed stellar energy) and cannot help 
to directly provide habitable conditions on its surface de- 
spite the low stellar flux. The enhanced heat flux (compared 
to that of the Earth) due to both the high planetary mass 
and the tidal heating might, however, have indirect effects. 
By inducing more volcanism and lasting longer over the 
planet's history, a higher heat flux favors the maintenance 
of the high CO2 level required on this planet. This may be 
important when taking into account that the system may 
be significantly older than the Earth. 
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4. Discussion 

4.1. Uncertainties in the location of the HZ limits 

Different uncertainties affect the determination of the 
boundaries of the HZ. Recently released data on H2O, such 
as the spectroscopy and equation of state at high pressure 
and the temperature and Raylcigh scattering cross-sections, 
could induce slight modifications on the planetary energy 
budget and should be included in radiative-convective mod- 
els. Other data are still lacking, including detailed param- 
eterization of the CO2-CO2 collision-induced absorption, 
which is an important factor for the greenhouse warming 
close to the outer limit of the HZ. 

For H 2 0- and CCVrich atmospheres, respectively at 
the inner and outer edge of the HZ, the sensitivity of the 
planetary albedo on the stellar effective temperature, from 
which Eqs. ^ and © are derived, was calculated assum- 
ing a blackbody spectral distribution of the stellar flux. In 
reality, TiO and H2O absorption makes the spectrum of 
M-type stars quite different from a black-body. The calcu- 
lation should be done for realistic spectra. For Gl 581, the 
expected effect would be to slightly move the HZ bound- 
aries away from the star, because TiO bands mostly block 
visible light and shift the emission towards the IR. The 
albedo sensitivity to T e g should also be studied for cloudy 
atmospheres of planets. Since clouds are a crucial factor in 
the location of the HZ boundaries, future simulations will 
have to include the effect of the spectral type of the star, 
in particular for CO2 clouds whose properties are very sen- 
sitive to the wavelength of the incident flux. 

An important step forward, especially in predicting the 
distribution of clouds, will be the use of 3-D global climate 
models (GC Ms), as was done to study syn chronously rotat- 
ing planets (|Joshi et al.lll997| ;[joshi 2003| ) or to add ress the 
runaway greenhouse effect ( Ishiwatari et alj |2002| ) . These 
models will have to include realistic microphysical processes 
for the formation, growth, and destruction of droplets and 
icy particles. To reduce the computing time, the treatment 
of radiative transfer in GCMs has to be greatly simplified. 
This is currently the main limitation to the application of 
GCMs to exoplancts. At the moment, the detailed line- by- 
line modeling required to accurately describe the transfer 
of radiation in planetary atmospheres is restricted to 1- 
D simulations. The development of 3-D GCMs, including 
robust treatment of radiative transfer, cloud physics, and 
photochemistry, represents the future direction of theoret- 
ical studies of habitability. 

4.2. Geophysical influences 

In Von Bloh et al. (|2007l . hereafter VB07), the estimate 
of the surface temperature is done by solving a system of 
coupled equations, including many geophysical processes 
affecting the CO2 atmospheric level (such as weathering 
and outgassing rates) and a radiation balance equation (Eq. 
4 of their paper), written as follows: 
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[1 -a(r surf ,P C o a )] =4/r(T s , 
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(6) 



where a is the Bond albedo, /r the outgoing thermal emis- 
sion and d the orbital distance. As shown in the present 
study, calculating a and Ir requires detailed radiative- 
convective modeling and, besides their dependence on T sur f 



and Pqo 2 > a an d Ir are highly sensitive to the spectral type 
of the star, the cloud properties and cover, and the abun- 
dance of other radiativcly active atmospheric compounds. 
Solving Eq. ([6]) numerically requires adopting a simple 
mathematical form for a(T surf , Pqo 2 ) an d ^R^surf, -FboJ, 
for instance a fit of climate model results. This mathe- 
matical form and the climate models from which it de- 
rives are not given in VB07, but most likely they use the 
same method as t he previous works of the group (e.g., 
iFranck et al1l2002t l. This is based on an interpolation of 
radiative-convective results obtained fo r spec ific cases by 
Kasting (1993) and Caldeira & Kasting (jl992t ) for a cloud- 
free atmosphere. This is consistent with the HZ inner limit 
shown in VB07 (figure 3b) corresponding (at t = Gyr) to 
a radiative-convective calculation with T sur f = 100°C and 
no clouds. There are several points we would like to discuss 
here: 

i) A parameterization derived from Kasting et al. (1993) 
can only be applied close to the conditions where the 
full radiative-convective modeling is available, which is 
close to the inner edge of the HZ. The radiative HZ 
would include Gl 581d only if the effect of CO2 clouds 
or another greenhouse gas is included. As there are only 
very few of these calculations at very low stellar irradi- 
ation, it would be very useful to know how a and Jr are 
calculated in VB07, in particular close to the outer edge 
of the HZ. 

ii) As shown in the present study, clouds have a major but 
uncertain effect on the location of both edges of the 
HZ. The inner edge, for instance, can be located any- 
where between 0.05 AU (runaway greenhouse with 100% 
cover of highly reflective clouds) and 0.11 AU (100°C, 
no clouds) from the star Gl 581. Therefore, the possible 
effects of clouds have to be compared with the geophys- 
ical influences included in VB07. 

Hi) Although the stellar luminosity is fixed in VB07 (which 
is a good approximation for M stars), one can see in 
their Fig. 3 that the inner boundary of the HZ (fixed 
at T sur f = 100°C) moves outward as the planet evolves. 
Unfortunately, this effect is not explained in the paper, 
so we do not know what geophysical process is involved. 

The present study and VB07 present two approaches 
that should complement each other. On one hand, we pro- 
vide atmospheric models developed to compute the radia- 
tive transfer and atmospheric structure in a self-consistent 
way. On the other hand, VB07 use a box model including 
geophysical processes relevant for the maintenance of habit- 
able conditions. Each of these processes has to be included 
in the VB07 model as a simple mathematical equation (e.g., 
Eq. [6]) to allow the numerical calculation of all the vari- 
ables through a system of coupled equations. Although we 
fully recognize the potential usefulness of this approach, we 
would certainly appreciate an investigation of the sensitiv- 
ity of the VB07 model to all the parameters involved and 
some estimation of the uncertainties affecting the model 
predictions. 

The most noteworthy effect shown in VB07 and previ- 
ous works from this group is the narrowing of the actual HZ 
compared to the radiative HZ presented for instance in our 
present study. The habitability of a planet can indeed be 
frustrated by the lack of outgassed CO2, as it may have oc- 
curred at some point in Mars' history. The age at which the 
decrease of internal heat and outgassing makes a planet no 
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longer habitable obviously depends on the planetary mass 
and composition. For a given stellar type, age, and plane- 
tary mass, the VB07 model computes the location of the 
frontier between the inner region of the HZ that can ac- 
tually be searched for habitable planets and the outer HZ 
region that can host only the dead remains of formerly hab- 
itable planets. However, to define this frontier accurately re- 
quires a robust model for each individual geophysical pro- 
cess, as well as strong constraints on the possible variety 
of planetary compositions; otherwise, we may rule out as 
"uninhabitable" circumstellar regions that could be prime 
targets in the search for habitable worlds. 
Two recent publications reach opposite conclusions on 
the onset of plate tectonics on larg e terrestrial planets. 
According to Valencia et al. (|2007af ). plate tectoni cs be- 
come more likely as planets get bigger. O'Neill et al. (|2007f ) 
come to the opposite conclusion that the higher gravity of 
big planets tends to prevent the formation of plates, pro- 
ducing a thick and unique crust. This debate illustrates 
the present difficulties in extrapolating Earth's geophysical 
models to exoplanets. Research on these important topics 
must of course keep going, but it is probably safer at this 
point to keep the broadest and less specific definition of 
the HZ boundaries from Kasting et al. (1993 J3 an d to allow 
future observations to determine whether a planet found 
within these boundaries actually is habitable. 

5. Conclusion: Are Gl 581c and Gl 581d good 
targets for Darwin/TPF? 

According to our present knowledge, based on available 
models of planetary atmospheres, and assuming that the 
actual masses of the planets are the minimum masses in- 
ferred from radial velocity measurements, Gl 581c is very 
unlikely to be habitable, while Gl 581d could potentially 
host surface liquid water, just as early Mars did. 

Because of the uncertainties in the precise location of 
the HZ boundaries, planets at the edge of what is thought 
to be the HZ are crucial targets for future observatories 
able to characterize their atmosphere. At the moment, our 
theory of habitability is only confirmed by the divergent 
fates of Venus and the Earth. We will have to confront our 
models with actual observations to better understand what 
makes a planet habitable. The current diversity of exoplan- 
ets (planets around pulsars, hot Jupitcrs, hot Neptunes, 
super-Earths, etc) has already taught us that Nature has 
a lot more imagination when building a variety of worlds 
than we expected from our former models inspired by the 
Solar System. 

It is obvious that the idealized model of a habitable 
planet atmosphere, where the two important constituents 
are CO2 and H2O, CO2 being controlled by the carbonate- 
silicate cycle, is likely to represent only a fraction of the 
diversity of terrestrial planets that exist at habitable dis- 
tances from their parent star. As an example, planets fully 
covered by an ocean may be common, either because they 
are richer in water than Earth or because the distribu- 

1 The precise location of the HZ boundaries based on this 
definition can of course be revised to account for new molecular 
data, additional greenhouse agents, 3D calculations, or better 
cloud models as done, for instance, by Forget and Pierrehumbert 
(1997), who estimated the possible effect of CO2 clouds in the 
location of the outer boundary. 



tion between surface and mantle water is different, or per- 
haps simply because, for a given composition, the mass-to- 
surface ratio and thus the water-to-surface ratio in crease s 
with the planetary mass, as noted by Lissauer (|l999h . 
Without emerged continents, it is not at all clear that the 
carbonate-silicate cycle could operate. The planets around 
Gl 581 can fall into this category since they are signifi- 
cantly more massive than the Earth (especially the >8 M s 
planet Gl 581d) and also because they may have started 
their formation in the outer and more water-rich region of 
the protoplanetary disk. 

Darwin/TPF-I and TPF-C could eventually reveal what 
the actual properties of the atmosphere of Gl 581c and 
Gl 581d are. From their thermal light curves we could infer 
if a thick atmosphere is making the climate more or less uni- 
form on both the day and night hemispheres of these plan - 
cts, despite a (nearly?) synchronized rotation (Selsis 2004). 
Visible and mid-IR water vapor bands could be searched in 
the atmosphere of Gl 581d to confirm its habitability. Mid- 
IR spectra of this planet could also reveal other greenhouse 
gases at work. Spectral observations of Gl 581c could po- 
tentially distinguish between a Venus-like atmosphere dom- 
inated by CO2 or an H2 0-rich atmosphere. The detection 
of O2 on this planet would generate a fascinating debate 
about its possible origin: as either a leftover of H2O photol- 
ysis and H escape or a biological release. There is certainly 
no doubt that Gl 581c and Gl 581d are prime targets for 
cxoplanct characterization missions. 
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